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The hadronic light-by-light scattering contribution to the muon g — 2 is the most troublesome 
component of its theoretical prediction; (1) it cannot be determined from the other measurable 
quantities, (2) the dimensional argument and the estimation based on hadronic models imply that 
the magnitude of this contribution may be comparable to the discrepancy between the standard 
model prediction and the experimental value. The direct approach to evaluate the hadronic light- 
by-light scattering contribution requires the evaluation of the correlation function of four hadronic 
electromagnetic currents, and the summation of it over two independent four-momenta of off-shell 
photons, which is far from the reach of direct lattice simulation. Here we propose an alternative 
method using combined (QCD + QED) lattice simulations to evaluate the hadronic light-by-light 
scattering contribution. 
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The muon anomalous magnetic moment (g — 2) has been measured very accurately at BNL [[I]] 

Aa p (EXP) =6x 1(T 10 . (1) 

Theoretical interest in the muon g — 2 stems from the fact that it can be used as a probe of unknown 
microscopic structures, such as supersymmetry with its breaking scale m§ ~ ^(100) GeV and 
large tan/3(> 10), and TeV scale gravity [gfl, each of which can give the additional contribution 
^(new physics) = <^(10 ~ 100) x 10~ 10 . We recall that both of these two structures are the most 
promising possibilities to solve the hierarchy problem, i.e., the hierarchy between the electroweak 
scale and the GUT or Planck scale. Thus, the currently available accuracy of Aa^ (EXP) can provide 
us with an opportunity to examine if the underlying structures ensuring such a hierarchy exist. 

Table [T] summarizes the current status of the muon g — 2 [Q], showing 2.7a deviation of 
tty (EXP) from the standard model prediction (SM) [[I], ^J. However, the uncertainty in (SM), 
which is attributed to the QCD contribution, must be reduced to elaborate such a discrepancy. 





x 10 10 


Aa^ x 10 10 


QED 


11 658 471.94 


0.14 


QCD 


695.4 


7.3 


Weak 


15.4 


0.22 


Theory 


11 659 182.7 


7.4 


Experiment 


11 659 208 


6 


a M (EXP)-a M (SM) 


25.3 


9.5 



Table 1: Comparison of a M (SM) with flu (EXP) 



Table ^ summarizes the QCD contribution to (SM) as well as its various components which 
are relevant in considering Aa jU (EXP) in Eq. (|J). The QCD contribution starts at ^(ofg m ) through 
the hadronic contribution to the vacuum polarization of the photon, which has been computed with 
lattice QCD [Q], and by using chiral perturbation theory and a vector meson model to fit the lattice 
data []5| . The improvement of the accuracy of this contribution either relies on the further precise 
measurement of o(e + e~ — ► hadrons)(y / i) [0], or on the improvement of the lattice calculation with 
a statistical uncertainty less than 5 % [Q] (but with less well-known systematic uncertainty). 





x 10 10 


Aa^ x 10 10 


hadronic vacuum polarization (^(cc^ m )) 


693.4 


6.4 


hadronic vacuum polarization (^(cc^ m )) 


-10.0 


0.6 


hadronic light-byTight {0{al m )) 


12.0 


3.5 


Total QCD 


695.4 


7.3 


Table 2: QCD contribution to the 


muon g — 2 


i 



Our target here is another type of QCD contribution, namely the hadronic light-by-light scat- 
tering (h-lbl) contribution to the muon g — 2 shown in Fig. [j], which also affects the physical 
interpretation of the muon g — 2 in viewing Table |2[ The diagram in Fig. [j] arises through the elas- 
tic scattering amplitude of two (off-shell) photons by QCD (the blob). To date, this contribution 
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could be estimated by purely theoretical calculation. So far, it has been calculated only based on 
the hadronic picture [Q, gj. Thus the first principle calculation based on lattice QCD is particularly 
desirable. 




fad 



elastic scattering amplitude 
of two photons by QCD 



Figure 1: hadronic light-by-light scattering contribution to the muon g — 2 

The diagram in Fig. [I] evokes the following naive approach; we calculate repeatedly the cor- 
relation function of four hadronic electromagnetic currents by lattice QCD with respect to two 
independent four-momenta l\, h of off-shell photons, and integrate it over l\, li- Such a task is too 
difficult to accomplish with use of supercomputers available in the foreseeable future. 

Here we propose a practical method to calculate the h-lbl contribution by using the lattice 
(QCD + QED) simulation; we compute 





QCD+quenched QED A \ / QCD+quenched QED g 

(2) 

quenched QED A 

amputate the external muon lines, and project the magnetic form factor, and divide by the factor 
3. In Eq. @ the red line denotes the free photon propagator D^ v (x, y) in the non-compact lat- 
tice QED solved in an appropriate gauge fixing condition. The black line denotes the full quark 
propagator Sf(x, y; U, u) for a given set of SU (3)c gauge configuration {U X)fl } and U (l) em gauge 
configuration {u x ^}, where the sum over relevant flavors / is implicitly assumed. The blue line 
represents the full muon propagator s(x, y; u). The average (, ) above means the one over the 
unquenched SU(3)c gauge configurations and/or the quenched t/(l) em gauge configurations 1 as 
specified by the subscript attached to it. Since two statistically independent averages over t/(l) em 
gauge configurations appear in the second term, they are distinguished by the labels A, B. 



'For the unquenched QCD plus quenched QED to respect the gauge invariance of QED, the electromagnetic charges 
of sea quarks are assumed to be zero. 
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To explain the mechanism underlying our method, let us look at the first term of Eq. 



perturbatively with respect to QED. Its magnetic components up to 



(a em ) consist of 



(«em) + ^(«em) 




QCD 



l+^(a cm ) 



(«em) 





QCD \ / QCDx<9'(a cm ) 



(3) 



The diagram in the first line gives the ^{cc^ m ) -contribution. The €?(a em ) -corrections to its muon 
part and to its quark part induce the ^(cc^ m ) -contributions shown in the first line and in the right 
diagram on the second line respectively. We recall that the QED gauge configurations in the first 
term of Eq. (^) are commonly shared by the quark part and the muon part. Hence, the photons 
can be exchanged between the two parts. As a consequence, the left diagram in the second line of 
Eq. (0) is induced at ^(o; e 3 m ), which takes the form of our target, Fig. |]. Alternatively, the quark and 
muon parts in the first and third diagrams in Eq. (^) are connected only by a single photon attached 
a priori. The second term in Eq. (Q) also contains those extra diagrams. Thus, by subtracting the 
second term from the first term, we may extract the h-lbl contribution. 

The quantities evaluated in our method (|J) are constructed from two currents for both terms, 
which are surely less noisy than the case of four currents encountered in the naive approach. Amaz- 
ingly, the only difference between the first and second terms of Eq. (||) is the ways averaging over 
the U(\) em gauge configurations. The h-lbl contribution should thus emerge as such a subtle dif- 
ference in averaging procedure. 

For our method to work efficiently, it is important to implement the cancellation of ^{cc^ m )- 
term in Eq. (|3j> without too much statistics. Sharing the same ensemble of U (l) em gauge configu- 
rations between the first term and the QED A part of the second term is one such method since the 
&{(xl m ) -correction to the muon vertex in both terms will be highly correlated. Furthermore, we 
organize Eq. (§) in the manner; 

H^(t F ,x c , U) 

= T ( ( LL 6 '^ u *) eYp^e^' s(x, xu u A ) 

\ W X F Xl 

xD pX (x, y) 

x HQf { tT {Yn s f( x ay> u > u A)ejxSf{y, x c ; U, u A )) 

f 



{tr (YpS f (x c ,y;U,UB)eYpSf(y,x c ;U,u B ))) )) ) , (4) 

7 ' "A I U 
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where the local form for the electromagnetic currents is assumed to avoid the complicate expres- 
sion, and the wave function renormalization constant of the muon is multiplied to remove the 
UV singularities in the vertex function of the muon. The ordering of the averages in Eq. @ is or- 
ganized in such a way that the cancellation of <^(ag m )-"components" can be realized as accurately 
as possible at the level of each configuration (U, ua). The nested sum for x, y could be efficiently 
carried out by working in the momentum space, where the Fourier transform of D(x, y) is known 
exactly. 

The observation done around Eq. (^) relies on the QED perturbation, while the simulation 
treats QED nonperturbatively. Thus, for our method to work, it is necessary that (a) the simulation 
can be done at the small coupling constant a em , possibly at its physical value a em = ^j, and that 
(b) the results for the weak coupling allow perturbative interpretation. As for (a), the quenched non- 
compact QED allows us to get the uncorrected C/(l) em gauge configurations efficiently As for 
(b), the lattice (QCD + QED) study of the electromagnetic splitting of the pion masses demonstrates 
that the splitting is perturbative JTc|], confidence that our method for the h-lbl contribution will 
work. 

So far we have focused on the diagram consisting of a single quark loop with electromagnetic 
currents inserted. In general, Fig. [j] contains two, three or four quark loops, where each quark 
loop contains at least one current. In the quenched approximation to QED, they must be computed 
separately. To illustrate implication of multi-quark loop diagrams, let us consider the charged 
pseudoscalar meson contribution [^] but with all electromagnetic currents put on its valence quarks, 
for simplicity. Three types of two quark loop diagrams provide the meson diagrams where the 
currents lie on both of two valence quarks, while the single quark loop diagram can provide the 
diagram where all the currents lie on one of two valence quarks. As is well-known, the calculation 
of multi-quark loop contributions (disconnected diagrams) is difficult in lattice calculations, and 
will be surely harder than the single quark loop diagram. Actually, the multi-quark loop diagrams 
vanish in the exact flavor SU (3) limit, except for a two quark loop diagram where four currents 
distribute over two quark loops evenly. It is an important issue to examine if they are subdominant 
and safely neglected. 

We thank RIKEN, Brookhaven National Laboratory and the U. S. Department of Energy for 
providing the facilities and hospitality where this work was done. M. H. 's work is supported in 
part by Grant-in- Aid for Scientific Research (15740173) in Japan. 
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